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Abstract. - We have investigated the collective modes of unconventional charge and spin 
density waves (UCDW, USDW) in quasi-one dimensional systems in random phase approxima- 
tion. The density correlator regains its normal state form due to the phase degree of freedom 
of the condensate. The possible effect of impurities is also discussed. From this, the current- 
current correlation function is evaluated through charge conservation. The spin susceptibility of 
USDW remains anisotropic an spite of the lack of any periodic modulation of the spin density. 
In UCDW, the spin response gets weaker in all three directions as the temperature is lowered. 



Unconventional density waves were intensively studied over the past few years. The so 
called d-density wave, which is UCDW with a gap of d x 2_ y 2 symmetry, has attracted much 
attention Q. This model successfully described many aspects of the pseudogap phase of 
high temperature cuprate superconductors. Also d-wave SDW can account for the anoma- 
lously small magnetic ordering in URu2Si2 |§, §]■ In quasi-one dimensional systems, our 
UCDW model describes successfully the behaviour of the low temperature phase of cv-(BEDT- 
TTF) 2 KHg(SCN) 4 [|[|,|], including the angular dependent magnetoresistance jjj. The re- 
sponse of (TMTSF) 2 PF 6 for T < T* ~ 4K can be described if the creation of USDW is 
assumed on top of the existing SDW § . 

The effect of collective modes is very important to study the dynamics of a system. The 
pole structure of the susceptibility is of prime importance J9|,[l^, [H]: at a given combination 
of the wavevector and frequency, the response function can be divergent. This means, that 
even without external field, excitations occur in the system whose dispersion is determined 
by the poles of the susceptibility. One way of identifying the correct order parameter among 
possible candidates is to study the unique collective modes supported by the ground state of 
a given symmetry |p^,p^,p^,[T^,[l^], as it was done in unconventional superconductors. 

Allowing for a low temperature phase transition, the Green's function of the model changes 
and new type of anomalous Green's functions can enter (for example the Gor'kov type 
F(ioj n ,k) function belonging to the (a J ^a^ k J) pair correlation in superconductivity). Con- 
sequently new types of interaction channels might open due to the novel type of nonvanishing 
expectation value, which is called the fluctuation of the order parameter As a result, the 
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simple RPA series is modified, and a number of coupled geometrical series are to be summed 



up. In the case of a conventional DW, the collective modes are well-known | IT , 18| , |lS| , [20|,|2l| , 
Generally one can conjecture whether the effect of RPA is or is not to re-establish the original 
pole structure of the metallic state. This can be achieved due to the degree of freedom of the 
phase of the density wave or the direction of the spin polarization of SDW. If we face with 
a magnetic phase transition (with a given spin orientation, i.e. SDW), the transverse (to the 
preferred direction) spin susceptibility and the density correlation function are not relevant 
from the transition's point of view, hence they will be restored. The former does due to the 
low energy (gapless) excitations of the spins in the perpendicular directions while the latter 
due to the free phase. Similarly the density correlator would regain its original form after 
RPA in a CDW due to the unrestricted phase of the density wave, but through the phonon 
mediated interaction the collective contribution is modified because of the mass enhancement. 
The quantities which are expected to change are the longitudinal spin susceptibility in SDW 
and all the spin susceptibilities in CDW because there are no low energy excitations based on 
the degrees of freedom of the system. 

From now on, we introduce the formalism necessary to investigate the collective modes in 
RPA and evaluate possible excitations of the main quantities. 

Formalism for RPA. - To start with, it is useful to introduce the spinor, which covers 
the whole momentum-spin space: 



*(k,T) 



/ a k,T(T) \ 

fl k-Q,T 



(1) 



where a^ ^ is the annihilation operator of an electron of momentum k and spin cr, Q is the 
best nesting vector. From this the Green's function of USDW is obtained as 

n( , • \ iuj n + t(k)p 3 + A(k)picr 3 

for UCDW (73 has to be replaced with 1, where <Ji (i = 1, 2, 3) are the Pauli matrices acting 
on spin space. With this, the interaction responsible for the UDW formation is determined 
from the gap equation of Ref. |2^] and assuming A(k) = Asin(bk y ) without loss of generality, 
it is given by 

N ~ 

— V(k, k', q, a, a) = (2 J y sm{bk y ) sin(6(fc^ - q y ))- 

~2F y sin(Wcy) sm(bk' y )) + 5 a ^ (J y - V y ) sin(bk y ) sin(b(k' y - q y )), (3) 

Of course the interaction is able to support Acos(6fc y ) type of gap, but we neglected the 
terms favouring this specific wavevector dependence: they are irrelevant with respect to RPA 
because they can only renormalizc the coefficients of the susceptibilities but are unable to drive 
the system into the desired ground state, namely with sinusoidal gap. This approximation 
was used in the first step of the calculation, namely in the gap equation p3f , when the specific 
wavevector dependence of the gap was chosen, because the matrix elements favouring this 
ordering are assumed to be the strongest. All the following calculations apply also to a 
cosinusoidal gap. For simplicity we shall limit our analysis to q = (q x ,0,0) (i.e. wavevector 
pointing in the quasi-one dimensional direction). 
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Density correlator and complex conductivity. For the density-density correlator we 
obtain for UDW 

{[n,n)) = ([n,n]) + ^{[n,A i ]) ([A i ,n]), (4) 
([Ai,n]) = ([Ai,n}) + ^-([A i; A i ]) ([A J , n]), (5) 

where (. . . )o means the thermal average without interaction between fluctuations. Here i = c, 
s for UCDW and USDW, respectively. A s — P2C3 sm(bk y ) 1 A c = p2sm(bk y ) and both P c 
and P s are positive. The detailed form of the interaction (P) can be found in Ref. p3j. The 
retarded products ([n,n])o, etc. are evaluated within the standard method and we find 

([n, n]) = 2.9(0)^-^ (l - 4A 2 F) , (6a) 
([n,Ai]) = -i4g(0)$AF, (6b) 
([A u Ai]) = 2.g(0) (-^ + (c 2 - e)F^j , (6c) 

where g(0) is the density of states at the Fermi energy in the normal state per spin, £ = vfQx, 
([n, n])o(q, uj) is obtained, for example, after analytic continuation from 

([n,n]} (q,iw„) = ^Tr(G(k,iw n )G(k - q, «<;„_,,)) (7) 



and 



k , n 



00 2tt 



- \; L . - - _j .tin i 1; 



^2 ,2 \ 2 ^ 7712 / /-2 1 , ,2 \ 1 yi A 2 _;__/'.,> 2 ^-2 



TV = (£ 2 - u 2 ) 2 - 4E 2 (e + u 2 ) + 4A 2 sm(y) 2 e, (9) 
D = N 2 — QAE 2 uo 2 ^ 2 {E 2 - A 2 sin(y) 2 ), (10) 

is the F function which also appears in the correlation functions of conventional DW with 
constant gap pcj j. The extra sin(y) 2 factor in the numerator of the density of states like term 
comes either from the gap or from the interaction. The remaining angle integral (which is the 
Fermi surface average) can be performed and after straightforward manipulation it yields to 



dE tanh — 
2 




(11) 
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where K(z) and Ii{n,z) are the complete elliptic integral of the first and third kind |24j| . 
respectively and 

N X = (-E 2 - e + (E - u) 2 )(e - (2E - lo) 2 ), (12a) 

Di = (E 2 + £ 2 — (E — uj) 2 ) 2 — 4£ 2 E 2 , (12b) 

N 2 = (~E 2 -£ 2 + (E + Lu) 2 )(f - (2E + u;) 2 ) : (12c) 

D 2 = (E 2 + £ 2 - (E + uj) 2 ) 2 -\i 2 E 2 . (12d) 

Putting these expressions together we obtain 

([n,n])(o;, g ) = 2^(0)^—2, (13) 

which is the same as in the normal state. It is worth mentioning that in conventional SDW 
the coefficient of £ 2 in the denominator is 1 + Ug(o) U is the on-site Coulomb repulsion. 
The strength of interaction is missing here due to its zero average over the Fermi surface 
(See Eq. (|)). In the absence of the DW pinning, all the modifications due to the change 
of quasiparticle spectrum are exactly canceled from the contribution of the fluctuation of 
the order parameter. The coupled RPA equations consist of two parts: the single particle 
contribution (one bubble contribution) which arises from the thermal excitation across the 
gap and possibly from a remnant portion of the Fermi surface, and the collective part from 
the motion of the density wave as a whole (due to the freedom of the phase) [f25[ . In the case 
of ([n,n]}, these two cancel each other to give back the normal state form. 

In real systems, impurities are always present. Henceforth the easiest way to incorporate 
the effect of impurities is to modify Eq. (^kj) as in Ref |2^] as 

([A*, Ai]) = 2g(0) (-^r - (e +lo 2 - c 2 )f) , (14) 

where lo v is the pinning frequency. This is the zeroth order phason propagator, and the 
inclusion of lu p ensures the presence of gapped phason mode. In general, impurities act 
differently on UDW [M than on conventional DW, but extended impurities are able to pin 
the phase of UDW, inducing finite pinning frequency Q. As a result, the density correlator 
reads as 

([n 7 n})(cj,q) = 2g(0K 2 f ~ ^t^fV" ^ ( 15 ) 

for both USDW and UCDW, / = 4A 2 i^. The zero sound dispersion (cj 2 = ^ 2 ) cancels out 
here contrary to the conventional case p(J , because a detailed study of the / function reveals 
that / = 1 + const(£ 2 — ui 2 )/A 2 in the limit of u> tends to £. The pole of the propagator 
describes the pinned dynamics of the USDW and UCDW condensate as uj 2 — lo 2 + £ 2 . Due to 
pinning, the condensate does not move below the threshold frequency uj p . The sound velocity 
in the presence of the electron-phonon coupling (g e _ p h) is given by 



C = C o Jl-gl_ ph lim([n,n])(0 )q ), (16) 

where Cq is the sound velocity without the electron-phonon coupling. First, in the absence of 
the DW pinning the sound velocity is temperature independent and reads as 

C = CqVT - A (17) 
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with A = 2g1_ ph g(0). Second, in the presence of pinning we obtain for UDW 

C = Cq\/1 — A(l — p s ), (18) 

where p s is the static condensate density [ p3[ , p s — lim 9 ^o lim w _,o /■ Since p s increases as 
the temperature decreases, we predict increasing sound velocity by lowering the temperature 
below its critical value T c for both UDW system. 

From the density correlator the complex conductivity in the chain direction is obtained 
using charge conservation as 

*( U ,q) = 2g(0) ie ^ w —^ w —^—^ . (19) 

From this the optical conductivity of UDW reads as Rea(oj) — 2g(0)e 2 irS(cu) in the pure case, 
which is essentially the same formula as in the normal state. Note that this formula applies 
for UCDW because there is no mass renormalization associated with lattice distortion. If 
conventional CDW is deduced from the attractive Hubbard model [^5[[ljJ, its conductivity 
coincides with those of a conventional SDW (using a phonon mediated attractive interaction, 
the mass gets renormalized in CDW). In the q — > limit the conductivity is the same as in 
conventional SDW, only the / function should be replaced by its version corresponding to 
UDW: 

*(«,, 0) = 2. 9 (0)- (f - 1 - -P^) (20) 

IUJ \ UJ Z — LO* J 

with /o = 4A 2 F g= o- It approaches the one in the normal state for u> ^> ui p for USDW and 
UCDW, while this statement is not true in CDW because m*/m ^> 1. The peak at zero 
frequency in the optical conductivity moves to the pinning frequency, suggesting that the 
primitive mimicking of pinning was successful. 

The conductivity for electric field applied in the z direction remains unchanged because the 
interaction is unable to dress the single bubble contribution due to the wavevector dependence 
of the velocity. In the y direction where the gap is developed, the following equation is to be 
solved: 

p 

{\jy,3y]) = {Uv,3y])o ~ J^(\Jy,3y])o(\jy,3y]), (21) 

where Pj = —P c , the matrix element responsible for the UCDW instability, and the current 
correlator without interaction between fluctuations is obtained as 

([j y J y ])o(u,q) = 4e 2 vlg(0)^—2(e -^7), (22) 

which gives the paramagnetic part of the total conductivity. Here the explicit wavevector 
dependence of the gap plays an important role. If the system possess a gap with cosine, the 
interaction is unable to dress the one bubble diagram, hence it gives the total paramagnetic 
part. The RPA equation is true only for a gap ~ sm(bk y ). Adding the diamagnetic term to 
it, the complex conductivity reads as 



a vyO 2 — V^wO _ a ri) 

7yV = 1 9(0)Pj (*yyO ~ 



(23) 
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where o~ n — —2g{0)e 2 v 2 /iuj is the normal state conductivity, cr yy o is the conductivity without 
RPA. Only Fermi liquid type corrections are present. Hence the excitation spectrum in the 
y direction is given by lo 2 = (1 + g(0)Pj(l — /o))£ 2 - This is very similar to the zero sound 
dispersion in the x direction in conventional SDW in the presence of pinning [ p0[ . The sign 
of Pj is negative for UCDW but for USDW it can be positive as well. These calculations 
verify our earlier assumption in Ref. p3| , that the optical conductivity in the perpendicular 
direction is given by the quasiparticle contribution, only Fermi liquid renormalization may 
occur. 

Spin susceptibilities. As to the spin susceptibility, the transverse one in USDW is 
expected to return to its normal state form similarly to the conventional case and to exhibit 
the trivial Goldstone mode. The RPA equations read as 

P 

([f7l,CTl]) = ([(7l,(7l])o + —([(7 1 ,p 1 (T 1 })o([p 1 a 1 ,(T 1 }), (24) 

p 

([piCTi,cri]) = ([piCTi,cri]) + —([pia 1 ,pia 1 ]) {[p 1 c7 1 ,a 1 }), (25) 



where 



[a^a^^^Lie-AAWF), (26) 
([a 1 ,p 1 a 1 ])o = 4g(0)iAujF, (27) 



[pi<Ti,Pi<Ti})o = 2g(0) { ^5 + (cj 2 - e)F ) . (28) 



Putting these together we obtain 

([oi,(7i]) = 2 5 (0)^L. (29) 

In the static limit, this expression reduces to the well known Pauli susceptibility. The lon- 
gitudinal susceptibility of USDW shows only Fermi liquid type renormalization, hence it is 
given mainly by the one bubble contribution, which coincids with the expression found for 
([n, n]) (q, iio v ) in Eq. ([?]). For UCDW the spin response is described with this expression in 
either direction. Although there is no obvious magnetic long range order in USDW ]23| ], yet 
it retains the spin anisotropy of SDW: in the static, long wavelength limit the longitudinal 
susceptibility reads as ([0-3, er 3 ]) = 1 — p s , which vanishes as T goes to zero. The pole belong- 
ing to the usual spin wave dispersion to 2 = £ 2 cancels out similarly to the case of the density 
correlator in the presence of pinning. 

Conclusion. We have studied the relevant correlation functions with random phase 
approximation. Again as in superfluid 3 He ||[l(J , the quantity whose correlator is investigated, 
often couples to the fluctuation of the UDW order parameter and it is necessary to handle 
these fluctuations on equal footing. Due to the unrestricted phase of the density wave, the 
density correlator regains its simple normal state form as it does in conventional spin density 
waves. In UCDW, the phason mass remains unrenormalized because the interaction from 
which this phase is originated, is electron-electron interaction and there are no retardation 
effects associated with electron-phonon interaction 22 , 19|. Based on a very simple concept 



the pinning of UDW is incorporated into the theory and the sound velocity C is examined in 
the presence and absence of the pinning. Without pinning C does not change compared to the 
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normal state value. In a pinned UDW the sound velocity increases upon entrance in UDW. 
Additionally the conductivity in the chain direction is obtained from the density correlator 
through charge conservation. It does neither change compared to the one in the normal state: 
in the pure system Reer(aj) <~ S(lo), which inspired the first pioneers in this field to identify 
the origin of the mechanism leading to superconductivity as CDW formation. 

By evaluating the current-current correlation function in the perpendicular directions, our 
previous assumption in Rcf. p3[ is verified: collective contributions do not show up and the 
quasiparticle contribution is enough to describe the electromagnetic response of the UDW. 

The spin susceptibility of USDW shows antiferromagnetic like anisotropy |^7],|2£| in spite 
of the lack of magnetic ordering, while in UCDW the magnetic response is influenced mainly 
by the quasiparticle contribution, leading to the freezing out of the static susceptibilities at 
low temperatures in all directions, as we predicted in Ref. p3f| , 

* * * 
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